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Atmosphere  dynamics 

Chapter- 1 


EQUATION  OF  MOTION 

1.1.  Inertial  or  absolute  coordinate  system: 

Inertial  coordinate  system  is  one,  which  holds  good  the  Newton’s  law  of 
acceleration.  The  Newton’s  law  states  that  the  applied  force  is  equal  to  the 
product  of  the  mass  and  the  acceleration  produced.  A Cartesian  coordinate 
system,  which  is  fixed  with  respect  to  the  fixed  stars,  is  very  nearly  an  inertial 
frame.  However,  the  one  that  moves  with  some  acceleration  is  called  the  ‘non 
inertial  coordinate  system’  in  the  space. 

1.2.  Rotating  Coordinate  System: 

Newton’s  law  depends  entirely  upon  the  nature  of  coordinate  system  relative 
to  the  surface  on  which  it  is  acting.  So  on  the  surface  of  the  earth,  this 
coordinate  system  also  rotates  with  the  earth  about  its  axis  of  rotation.  As  a 
result  the  departures  of  acceleration  under  such  a rotating  coordinate  system 
have  to  be  worked  out.  The  coordinates  on  the  rotating  system  of  earth  taken 
are  positive  ‘x’  as  latitude  towards  east,  positive  ‘y’  along  the  meridian 
towards  the  North  Pole  and  ‘z’  as  the  height  into  the  atmosphere  as  shown  in 
Fig.l. 
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Fig.l.  Coordinates  on  the  spherical  earth 

1.3.  Equations  in  a rotating  Coordinate  System: 

Fet  us  consider  a coordinate  system  with  the  coordinates  x’,y’  and  z’  as 
shown  in  Fig.2.  Fet  the  coordinate  system  rotated  by  an  angle  Ut  and  the  new 
coordinates  are  x,y  and  z respectively.  From  the  rotational  coordinate  system 
we  can  write 
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From  triangle  ORS,  OR  = 


From  triangle  PQR,  QR  = OQ-OR  = X ; PR  = ^ A 

CosCit  cosClt 


x 

From  triangle  PQR,  SinCit  = ^ = C°sCLt  ^>X'=  XCo£lt  - TjSinflt 

cosQ  t 

Similarly  from  triangle  PTU,  PT  = X,  From  triangle  OUV,  TU  = 


_n_ 

cosQt 


v OT  = rj,OU  = - 


From  triangle  PTU,  TanTlt  = 


TU  _ cosQf 
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v[  = X Sin  Dt  + r|  Cos  Qt 

It  may  be  noted  that  the  letter  ‘ Y’  may  be  read  in  place  of  r|  then  the  we  get  the 
following  equations. 

X;  = X Cos  Qt  -YSinQt 

Y;  = X Sin  Ht  + Y Cos  fit  (1.1) 

7!  = Z 

Similarly  in  place  of  X,  Y,  and  Z if  forces  Fx,  Fy,  and  Fz  are  considered  we  can 
write 

Fx;  = Fx  Cos  Qt  - Fy  SinDt 

F/  = Fx  Sin  Qt  + Fy  Cos  Qt  (1.2) 

Fz;  = Fz 

If  we  differentiate  equation  1 . 1 twice  with  respect  to  time  we  get 


X'  = CosOt(X  X - 2yO  - XO2  - SinQt(y  + 2x0.-  yO 2) (1.3) 

According  to  Newton’s  law,  a = F/m  = X (1.4) 

So  equation  1.2  and  1.4  combining  we  write 
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since  X'  = Equations  1.3  and  1.5  can  be  equated.  After  equating  and 
comparing  the  coefficients  of  Cos  Dt  and  Sin  Dt  we  get  as  below: 


X-2YX-X02  = F*/ 
/m 

or  X = XQ2  + 21X  + F'/ 
/m 

Y + 2X0-Y02  =Fy 
/m 

Y = Y02  -2XO+Fy/ 

/m 

(1.6) 

Similarly  can  be  written 

Z-^  = 0 

m 


Equations  (1.6)  are  Newton’s  second  law  transformed  to  a rotating  coordinate 
system. 

The  first  of  the  terms  on  the  right  hand  side  (XO2,  YO2)  denote  the  centrifugal 
accelerations  pointing  radially  outward.  The  second  terms  on  the  right  hand  side 
(2  YX  ,2X0 ) represent  accelerations  owing  to  the  combined  effects  of  rotation  of 
the  coordinate  axes  and  motion  of  the  particle,  relative  to  the  rotating  coordinate 
system.  These  are  called  Coriolis  accelerations.  The  action  of  this  Coriolis 
acceleration  is: 

1 . It  acts  only  perpendicular  to  the  plane  of  the  moving  particle. 
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2.  It  acts  right  hand  side  in  the  N.H  and  left  hand  side  in  the  S.H  of  the 
moving  particle. 

3.  It  comes  into  being  only  for  a moving  particle  but  not  for  a fixed  particle. 

4.  It  does  not  change  the  magnitude  of  the  velocity  but  changes  only  the 
direction  of  the  moving  particle  as  it  acts  perpendicular  to  the  plane  of 
moving  particle. 

Thus  the  centrifugal  and  Coriolis  accelerations  are  extra  terms  obtained  due  to  non 
inertial  (non  Newtonian)  frame  of  reference. 

1.4.  Spherical  coordinate  system: 

We  shall  now  apply  these  results  on  a rotating  spherical  earth,  in  which  x is  tangent 
eastward,  y meridionally  northward  and  z radially  upward  into  the  atmosphere  as 
shown  in  Fig.3.  Let  ‘P’  be  a point  on  the  earth’s  surface  which  is  rotating  with  the 
angular  velocity  Q,  at  a distance  of  ‘acoscp’  from  the  axis  of  rotation  where  ‘a’  is  the 
radius  of  the  earth  and  ‘cp’  is  the  latitude.  From  fig.3,  Qx  = D Cos  90  = 0,  fly  = fl  Sin 
(90-cp)  = D Cos  cp,  Dz=  Q Sincp. 


Fig.3  a &b 

If  the  velocity  of  the  particle  is  V=  iu  +jv+kw  then  the  Coriolis  acceleration  is  curl  of 
the  vector  as  shown  in  fig.  3b.  The  direction  of  the  vector  D is  such  that  when  the 
fingers  of  the  right  hand  bend(curl)  about  the  axis  of  rotation  in  the  same  sense  as  the 
motion,  the  thumb  points  in  the  direction  of  H. 

In  fig.3b  let  ‘a’  is  the  position  vector  of  the  moving  particle’P’  as  measured  from  the 
origin  of  coordinates  (centre  of  the  earth).  So  the  vector  Dxa  has  the  magnitude  ‘D  a 
sincp’,  which  is  just  the  linear  speed  of  the  particle  (V).  Therefore  V = Dxa.  Similarly 
if  the  magnitude  of  velocity  (Y)  is  considered  as  the  position  vector,  then  the 
acceleration  is  Qxy.  Thus  we  got  the  Coriolis  accelerations  in  equation  1.3  as 

- 2XDand  - 2YQ. , where  XandY  are  the  velocities  ( — , — ).  Which  means  the 

dt  dt 

Coriolis  accelerations  on  the  surface  of  the  earth  can  be  written  as  -2Qx  y. 
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C = -2Dxy  = -2D 

i j 

0 cos^ 

k 

sin  (/) 

= -2D 

i 

Qx 

j 

Dv 

k 

D. 

u V 

w 

u 

V 

w 

where  u,v,w  are  the  x,y,z  components  of  the  velocity.  Expanding  the  determinant  and 
separating  the  different  components  we  get, 

Cx  = -2D  w coscp  + 2Dv  sincp 

Cy  = -2Dusin(p  (1.7) 

Cz  = 2D  u cos  cp 

The  Newton’s  second  law  of  motion  on  a rotating  coordinate  system  becomes  from 
equation  1.6  as 
..  F 

X = 2Q.vSiruf>  - 2£lwCos(j) 

m 

..  F 

Y — -2Q.uSin</>  (1.8) 

m 

..  F 

Z — f^p  2 Q.uCos</> 
m 

1.5.  Other  forces  or  accelerations: 

1.5.1.  Gravitation  and  Gravity: 

The  term  gravitation  solely  means  the  attractive  force  between  bodies  which  is 
described  by  Newton’s  law  of  gravitation.  This  force  of  attraction  is  proportional  to 
the  product  of  the  masses  and  inversely  related  to  the  square  of  the  distance  between 
them. 

ga  = -(GM/r2)r 

Where  G=  6.658xl0'8,  M = 5.988xl027  grams  is  the  mass  of  the  earth,  r = the  position 
vector  of  the  unit  mass  of  the  air  molecule  measured  from  the  centre  of  the  earth. 

As  the  particle  on  the  surface  of  the  earth  is  rotating  along  with  the  earth,  it  has 
centrifugal  acceleration  and  the  net  effect  appears  as  acceleration  due  to  gravity  as 
shown  in  Fig.4. 
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Fig  .4  Gravitational  and  centrifugal  accelerations  on  the  earth 

g = ga  - ^x(Qxr) 

which  is  equal  to  acceleration  due  to  gravity  = gravitation  - centrifugal  acceleration. 
At  the  surface  of  the  earth  acceleration  due  to  gravity(g)  depends  on  the  latitude.  At 
poles  it  is  maximum  and  minimum  at  the  equator. 

1.6.  The  Pressure  gradient  force 

Consider  a rectangular  volume  having  the  sides  8 x , SyandSz  fixed  in  a coordinate 
system  relative  to  the  solid  earth  as  shown  in  Fig.5.  Then  the  force  in  the  x-direction 
due  to  the  atmospheric  pressure  is,  p SyS,  where  ‘p’  is  the  atmospheric  pressure  and 
Sy8,  is  the  area  of  the  face  it  is  acting. 


Fig.5.  Pressure  force  acting  on  a rectangular  volume  element 
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Similarly  the  pressure  force  acting  on  the  opposite  face  would  be  - p-\ Sx  \dydz 

K dx  ) 

The  net  force  acting  in  the  x direction  (i)  would  be 

(ipSySz)  + \-i(p  + — Sx)SySz.  \ = -i  — SxSySz, , pressure  force  per  unit  volume=  - 
[ dx  ) dx 

iop/ox  because  the  volume  of  the  cube  is  Sxdydz. 

Similarly  considering  in  other  directions  the  total  pressure  force  would  be 

-\i?P  + j?E+k^\ 

\ dx  dy  dz\ 

The  force  per  unit  mass  = ~ — \i  — + j — + k — \=  -aVp (1.9) 

p [ dx  dy  dz] 


u 
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Fig.6.  As  the  pressure  increases  top  to  bottom,  pressure  force  acts  in  the  opposite  direction 

The  minus  sign  indicates  that  as  the  pressure  increases  in  one  direction,  the  pressure 
force  acts  in  opposite  direction  as  shown  in  Fig.6.  Because  of  the  presence  of  gradient 
(del,V)  operator,  it  is  called  as  the  pressure  gradient  force. 
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1.7.  FINAL  EQUATION  OF  MOTION: 


Finally  including  all  the  forces,  the  vector  equation  of  motion  can  be  written  as  ( from 
equations  1 .7) : Table  2.11 


Total 

acceleration 

Coriolis 

acceleratio 

n 

Pressure 

gradient 

Acceleration 
due  to  gravity 

Friction 

X = 

d2x  du 
dt2  dt 

f.v- 

2Dwcoscp 

dp 

-a— 

8x 

0 

Fx/m 

y _d2y  _ dv 
dt 2 dt 

-f.u 

'“I 

0 

Fy/m 

2 _ d2z  _ dw 
dt2  dt 

+2Ducos(p 

dp 

-a— 

dz 

-g 

Fz/m 

Then  we  can  write  the  Navier-  stokes  equations  as: 


Du 

~Dt 

Dv 

~Dt 

Dw 

~Dt 


i dp  / d2u  d2u  d2u\  n 

— — + W t-7  + TT  + TT  +(20usin<£  — 2nu;cos<£) 
p ox  \ ox 2 ay2  dz2 ) 

1 dp  (Pv  d2v  d2i A ^ 

--Pfr+v(w+w  + w)-(-2Slusm,l,) 


1 dp 

— ~f  + v 

p dz 


/ d2w  d2w  d2w\ 

\dx2+~df+~d^) 


+ (2Hmcos<^) 


\ 


acceleration 
of  air  parcel 


pressure  viscous  force 
gradient 


gravity 


t 

Coriolis  force 


(Or) 


du  du  du  du 

— + U—  + V—  + W— 

dt  dx  dy  dz 


dv  dv  dv  dv 

~r  + «T-  + V—  + XV— 

dt  dx  dy  dz 


1 dp  / d2u  dzu  d2u\ 

p dx  "I"  \ dx2  dy2  dz2  ) 


+ (20u  sin  (f>  — 2Qwcos  (f> ) 

1 dp  / dzv  d2v  <92v\ 

p dy+  \ dx2  dy2  dz2  ) 


— (20m  sin  <t>) 


dw  dw  dw  dw  1 dp  ( <92Jn  <92iu  d2w\ 

i;+uj;+vty+wj;=-g-pTz+v{^+~3?+i?) 


+ (2  cos  <f>) 


This  set  of  equations  is  usually  called  the  Navier-Stokes  equations  for  the  conserva- 
tion of  momentum.  The  equations  are  named  after  Claude  Louis  Navier  (1785-1836) 
and  George  Gabriel  Stokes  (1819-1903),  who  both  contributed  to  the  development 
of  the  equations.  Navier  developed  the  form  of  these  equations  for  an  incompressible 
(that  is,  constant  density)  fluid  in  1821.  In  1 822,  he  published  a further  refinement 


1.8:  SCALING  THE  EQUATION  OF  MOTION: 

Scaling  is  useful  to  know  which  term  weighs  more  than  the  other  terms  in  the  equation 
of  motion  and  so  the  terms  that  weigh  more  can  be  retained  and  the  others  can  be 
neglected  or  ignored  for  different  calculations.  The  procedure  is  to  take  approximate 
rounded  off  values  (order  of  magnitude)  of  each  term  and  substitute  as  per  the 
dimensional  value  of  the  each  term  as  given  in  the  table  below: 

1.8.1.  Scaling  the  different  terms  of  the  equation  of  motion  in  the  case  of  tropical 
weather  systems: 


Item 

Order  of  magnitude 

1.  horizontal  velocity  u ,v 

10  m/s  = 103  cm/s 

2.  vertical  velocity 

1 cm/s 

3. horizontal  length  of  a 
synoptic  system, Ax, Ay 

1000  km,  =10*  cm 

Vertical  depth  of  the 
system,  Az 

10  km,  = 10b  cm 

5.  pressure  fluctuation,  Ap 

10  mb  = 104  dynes/cm2 

6.time  dimension 

105  seconds 

1.  Acceleration:  du/dt  = 103/105  = 10'2 

2.  Coriolis  terms: 

a)  2 Dv  sin  cp  = 2x7.292xl0'5xSin90°xv  = 10xl0'5xlxl03  = 10'1  So  2D.  v simp  is 

valid. 

b)  2Qw  cos  cp  = 2x7.292x10'5xCos(0)xw  = 10xl0'5xlxl  = 10'4. 
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Which  means  2Gw  Coscp  is  negligibly  smaller  than  2GvSin(p.  So  2QwCoscp 
can  be  neglected. 

c)  2QuSincp  = 2x7.292xl0'5xSin90xu  = 10xl0'5xlxl03  = 101,  valid. 

d)  2Du  Cos  cp  = 2x  7.292xl0'5xSin90xu  = 10xl0'5xlxl03  = 101 

3.  Gravity  term:  g = 9.8  = 1 0m/s2  = 1 03  cm. 

4.  Pressure  gradient  terms: 

adp/dx,  adp/dy  = (1/p)dp/dn  = 10 4/  (10'3  xlO8)  = 10'1 
adp/dz  = (l/p)dp/dz  = 104/  (10'3xl06)  = 101 

i)  Scaling  x equation:  (frictionless  flow): 

— = -a  — + 2Qvsin®-2Qwcos®  = 10'2  = 10'1 +10'1  - 10'4  The  last  term  is 
dt  dx 

negligibly  small  compared  to  the  other  two  terms  in  the  right  hand  side  and  so  can  be 
neglected. 

dii  dp  __  . , 

.’.  — = -a  — + 2Qv  sin  O 
dt  dx 


ii)  Scaling  y-equation: 


= -a  — ~ IQuSind)  — > 10'2  = 10'1  - 10'1,  both  terms  in  the  right  hand  side  are 

dy 


dt 


equal  and  so  all  are  important. 

iii)  Scaling  z-  equation: 

dw  dp 
x — 
dz 

gradient  and  gravity,  2Qu  coscp  is  negligibly  small. 


= -a  — + 2D.u  cos®  - g — * 10'5  = 101  + 10'1  - 103,  when  compared  to  pressure 
dt  dz 


dw  dp 

— = -a  — - 2 
dt  dz 

Thus  final  scaled  equation  of  motion  is 


du  dp  _ 

.-.  — = -a  — + 2Qv  sin  ® 
dt  dx 

— = -a  — - 2CiuSinO 
dt  dy 

dw  dp 

— = -a  — -g 
dt  dz 
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Chapter-2 


2.  HYDROSTATIC  EQUATION  & GEOPOTENTIAL: 

The  third  equation  of  motion  reduces  to  hydrostatic  equation  as  dw/dt  and  2Du  coscp 
are  neglected. 

Therefore,  -adp/dz  - g = 0, 

= pg  o adp  = gdz^>  J adp  = J gdz  =^>  J d<t>  = J gdz 

This  cp  is  called  the  geopotential  and  has  the  dimensions  L2T'2.  The  unit  of 
geopotential  is  dynamic  meter  which  is  numerically  2%  smaller  than  the  geometric 
meter  because  while  the  work  needed  to  lift  a unit  mass  through  one  geometric  metre 
is  about  0.98  dynamic  meters,  which  is  close  to  1.0  geopotential  meter.  The  value  gdz 
is  the  potential  energy  between  the  thicknesses  of  two  levels.  Constant  geopotential 
lines  are  called  equigeopotential  lines  which  are  also  level  surfaces.  To  move  along  an 
equigeopotential  line  no  extra  energy  is  required  as  the  acceleration  due  to  gravity  is 
constant  on  the  equipotential  lines.  The  equigeopotential  surfaces  on  the  earth  are 
oblate  spheroidal  like  that  of  the  earth  as  equatorial  radius  is  larger  than  the  polar 
radius. 

2.1.  THE  EQUATION  OF  CONTINUITY: 


An  important  restriction  upon  the  velocity  of  a fluid  may  be  obtained  by 
kinematic  methods  applied  to  the  law  of  conservation  of  mass. 

Consider  an  infinitesimal  rectangular  volume  fixed  in  space  with  sides  of 
lengths  5x,  5y,5z  in  a moving  fluid  having  velocities  u,v,w  respectively  as  shown 
in  Fig.7.  Consider  first  the  flow  parallel  to  the  x-axis.  As  ‘u’  is  x-  component  of 
velocity  and  p is  density  of  the  fluid  the  mass  flow  through  the  left  face  is  p u. 
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The  mass  flow  into  the  volume  per  unit  time  is  ( p u)5y5z 

The  mass  flow  out  of  the  volume  per  unit  time  on  the  right  face  is 

( p + —Sx)(u  + — 8x)dydz 
dx  dx 

So  the  net  flow  out  of  the  volume  in  the  x-direction  is 

(y oudydz)  -\(p  + — Sx)(u  + — <Sr)ldy<% 

L dx  dx  J 

[ d dp  du 

I dx  y dx  dx 


As  the  volume  of  the  rectangular  element  considered  is  infinitesimal,  5x  is 
very  small  and  so  higher  order  term  (second  term)  is  negligible. 

Therefore,  the  net  flow  in  the  x-direction  is  = - — ( pu ) 

\_dx  \ 

Taking  the  mass  flow  in  all  the  three  component  directions,  the  total  flow 

out  would  be  - — (pu)  + — (pv)  + — (pw) 

| dx  dy  dz 

where  V and  ‘w’  are  the  velocity  components  in  the  y and  z directions 
respectively. 

We  shall  compute  the  rate  at  which  mass  enters  the  box  as  a result  of  fluid 
motion  and,  since  mass  must  be  conserved,  we  shall  equate  this  to  the  rate  of 
change  with  time  of  the  mass  contained  within  the  cube. 
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In  other  words,  if  mass  is  neither  created  nor  destroyed  and  if  the  sum  is  not 
zero,  it  must  be  balanced  by  a change  of  density  (rate  of  change  of  density  — ) 


-^-  = -f{pU)+-f(pV)  + ^-(pW)  = V.pV — (2.1) 

at  ox  dy  oz 

, ./,'  JVr  . ,t  . , >v  l 

p dt  lex  dY  ez\ 


dt 


This  is  called  the  equation  of  continuity.  The  first  term  is  the  fractional  rate 
of  change  of  density  and  the  second  term  is  the  fractional  rate  of  change  of 
volume. 

If  the  fluid  is  incompressible,  the  density  is  uniform  (constant),  then  ~ = 0 


p dt 


= 0 Then  the  equation  of  continuity  is 


^A  + ^ = v.v  = o... 

dx  dy  d z 


..  (2.3) 


where  V = iu  + jv  + kw 

Equation  (2.3)  is  frequently  considered  convenient  as  it  is  nothing  but 
horizontal  three  dimensional  divergence. 


Chapter  -3 

Circulation  and  vorticity 


3.1.  CIRCULATION: 

Let  the  fluid  is  flowing  along  a path  ‘ s’,  a portion  of  it  is  say  ‘ds’,  with  the 
components  dx  and  dy,  1 = cos  (3  = dx/ds,  cos  0 = dy/ds  = m;  where  1 and  m are  the 
directional  cosines  of  ‘s’.  Therefore  dx  = 1 ds  and  dy  = m ds  as  shown  in  Fig.8. 
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Fig-8 


We  know  from  potential  line  (orthogonal  stream  line)  equation  as  udx  + vdy  =0 


Substituting  in  above  for  ‘u’  and  V as 


dx  _ dx  ds  _ i 
dt  ds  dt 

dy  dy  ds  _ 
dt  ds  dt 


Therefore,  (lc.l  + mc.m)  ds  = 0 

(l2  + m2)  c ds  = 0 , so  C.ds  = 0 ( since  l2  + m2  = 1) 

or  udx  + v dy  = cds  = 0 (3.5). 

s2  d 

If  this  equation  is  integrated  between  Si  and  S2  [ (u  — + v—)ds  = \ C.ds , 

• ds  ds  J 


which  is  the  flow  of  fluid  between  Si  and  S2.  It  has  the  dimensions  cm2/sec 
( I?T~^ ).  As  the  velocities  u and  v are  each  functions  of  x and  y in  the  fluid  to 
perform  the  integration  it  is  necessary  to  know  what  these  functions  are.  They  can 
be  specified  along  a given  line  of  flow  and  the  resulting  integral  along  this  line  is 
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called  a line  integral.  Similarly  a line  with  angular  turns  is  called  a closed  line 
integral. 

So  the  flow  around  a closed  path  or  the  line  integral  around  a closed  path  is 

called  the  circulation ‘C’ = | udx  + vdy (3.6) 

Thus  ‘C’  is  a measure  of  the  extent  to  which  the  fluid  exhibits  rotary  motion. 

3.2.  VORTICITY: 


To  obtain  a quantitative  understanding  of  the  rotating  fluids  with  a change  of 
co,  express  the  motion  in  Cartesian  coordinates  of  x and  y.  The  rotation  is 
considered  for  two  infinitesimal  lines  dx  and  dy  in  ‘x’  and  ‘y’  directions.  Let  ‘u’  and 
‘v’  are  the  linear  velocity  components  in  y and  x directions  respectively  as  shown 
in  Fig.9a.  As  the  angular  velocity  is  different  on  the  two  axes 


v = coi  x and  -u  = -C02y  or 


8u 

- — = o)2 

dy 


2 


The  average  of  the  two  angular  velocities  say  ‘co’  is  co  = (coi+co2)/2  = 

(3j> 


where  is  called  the  vorticity  which  is  twice  the  local  average  angular 
velocity. 


/K 

^ Jk- 


't'  V 


Fig.9a 

For  cyclonic  (anticlockwise  in  northern  hemisphere)  flows  £ is  positive  and 
for  anticyclonic  flows  (clockwise  in  N.H),  £ is  negative  as  shown  in  Fig.9b. 


15 


(A) 

^IrjLpSl^ 


Fig.9b 

3.3.  PLANETARY  VORTICITY: 

In  the  case  of  earth  rotation  this  is  considered  as  vertical  component.  We 
know  that  the  rotation  of  the  earth  about  its  axis  results  in  the  deflection  of  winds 
by  Coriolis  force.  These  deflections  are  explained  in  terms  of  the  poleward  decrease 
in  the  eastward  velocity  of  the  surface  of  the  earth.  In  addition  to  the  linear 
eastward  velocity,  the  surface  of  the  earth  has  an  angular  velocity.  Like  the  linear 
eastward  velocity,  this  angular  velocity  depends  on  latitude.  At  the  poles,  the 
angular  velocity,  of  the  surface  of  the  earth  is  simply  O (=  7.292  x 10 5 rad  sec1).  At 
lower  latitudes,  it  is  a proportion  of  O.  The  larger  the  angle  between  the  earth’s 
axis  of  rotation  and  the  local  vertical,  the  smaller  is  the  angular  velocity  of  the 
surface  of  the  earth  as  shown  in  Fig.  10.  At  the  equator,  where  a vertical  axis  is  at 
right  angles  to  the  axis  of  rotation  of  the  earth,  the  angular  velocity  of  the  surface  is 
zero.  Thus  the  angular  velocity  of  the  surface  of  the  earth  about  a vertical  axis  at 
latitude  O is  given  by  O sin  ®.  As  the  vorticity  is  twice  the  local  average  angular 
velocity  (as  given  by  equation  3.7),  the  planetary  vorticity  is  equal  to  2(0  sin  O)  (= 
f).  The  planetary  vorticity  at  the  poles  and  the  equator  respectively  are  ‘20’  and 
zero. 
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Fig.  10. 


3.4.  ABSOLUTE  VORTICITY: 


We  then  define  the  absolute  vorticity  as  the  sum  of  the  two  vortices  ,(  £ + f ). 
The  vorticity  ‘C’  is  measured  rotation  of  the  particle  relative  to  the  surface  of  the 
earth  and  so  it  is  called  relative  vorticity. 

Considering  the  frictionless  horizontal  equation  of  motion: 


du  dp  , 

— = -a  — + iv 

dt  dx 

dv  dp  _ 

— --a— — fu 

dt  dy 


if  these  are  cross  differentiated  and  subtracted  one  from  the 


other,  to  eliminate  pressure  gradient  terms,  we  get:  — (C  + /)  = ~(C  + /)  — + — 

■"  A*  1 dx  dy 


Since  we  know  the  equation  of  continuity  ^ ^ = 0 , ^-  (£"  + /)  = 0 
upto  the  column  of  thickness  ‘D’  in  the  atmosphere  called  gradient  wind 

level  is  considered,  then  — + ^ ) = 0 . The  term  ^ + f is  called  the  potential 

dt  D D 

vorticity. 
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3.5.  VORTICITY  THEOREM: 


While  circulation  is  a measure  of  tendency  of  fluid  rotation,  vorticity  is 
the  same  property  of  fluid  for  an  infinitely  small  area. 

To  obtain  the  vorticity  theorem,  the  horizontal  equations  of  motion  are 


du  du  du  du  r dp  _ 

h u — + v — + w — = jv  -a l -Fx.. 

dt  dx  dy  dz  dx 

dv  dv  dv  dv  , dp  _ 

— + u — + v hw — = -ju  - a — + Fv . 

dt  dx  dy  dz  dy 

. . „ dv  du 

and  we  know  C = 

dx  dy 


..3.14 


..3.16 


d 

dv 

du 

d 

'dv 

du 

_ d 

dv 

du 

dt 

dx 

dy_ 

' U dx 

dx 

dy  _ 

' V fy 

dx 

dy 

-ich 


Take  partial  derivative  with  respect  to  (w.r.t)  y to  equation  3.14  and 
differentiate  w.r.t  x to  equation  3.15  and  subtract  3.14  from  3.15.  and  re 
arrange,  then  we  get 

d_ 

V dz 

dz  dx  dz  dy 

dFy  QFX 
dx  dy 

df 


'dv 

du 

du 

dv 

du 

dv 

_dx~ 

dy  _ 

+ dx 

dx " 

dy  _ 

dy  _ 

y_ 


=./r^+^i_La+vai_a[0]-te^-§p^i+r^L-^i 

\_dx.  dy]  |_  dx  dy]  dx  dx  dy  J [_  dx  dy  J 

v—  + w— 
dy  dz 
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A 

a* 


rdvdw  du  dw~\  _ df  I ' dp  da  dp  da]  I”1 

| _dz  dx  dz  dy\  dy  |_dy  dx  dx  dy  \ |_ 

(3.18) 


dFx 

dy 


dx  dy  dt  dy 

It  can  be  further  written  as 


— (r)+  —(f)=  -D(f  + C)—  ^ ^ du  dw  dp  da  dp  da  dFy 
dt dtKJ  ’ y ' ldz~d^c~^z^y\ldy  + 


d />  | _ jytf  \ /-)  l dv  dw  du  dw\  \dp  da  dp  da]  \dFy  dFx  1 

dt  |_Sz  dx  dz  dyj  |_dy  dx  dx  dyj  \dx  dy  J 

(3.19) 

I II  III  IV 


This  is  called  vorticity  equation.  The  term  on  the  left  hand  side  ( r.h.s)  is 
called  the  rate  of  change  of  absolute  vorticity.  This  equation  says  that  the  absolute 
vorticity  of  a parcel  of  air  can  change  only  through  the  contribution  of  four  terms 
on  the  right  hand  side. 
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Show  that  the  equation  for  the  time  evolution  of  £is 

dt_ 

-(/+£) 

du  dv 

dwdv  dwdit^ 

dp  dp 

dp  dp] 

./I 

dy 

dt 

,dx  dy, 

Kdx  dz~  dy  dz, 

k dx  dy 

dy  dx) 

convergence  tilting  Baroclinic  term  “beta”-term 

...,3.19b 

3.6.  PHYSICAL  INTERPRETATION  OF  THE  TERMS  OF  VORTICITY  EQUATION: 

The  first  term  (I)  on  the  r.h.s  of  equation  (3.19)  is  called  the  divergence  (or 
convergence)  term.  Positive  divergence  increases  the  radius  of  rotation.  As  the 
term  is  negative,  it  decreases  the  magnitude  of  absolute  vorticity. 

The  second  term  (II)  of  equation  3.19  on  the  r.h.s  is  known  as  the  tripping 

dw 

(or  tilting)  term.  When  the  vertical  shear  ( — ) of  the  first  half  of  this  term 

dx 

decreases,  horizontal  shear  ( — ) increases  upward.  As  it  has  minus  sign,  its 

dz 

contribution  to  (f+  C ) is  positive.  Which  means  (f+  Q will  increase  with  time.  A 
corresponding  analogy  can  easily  be  made  for  the  second  half  of  this  term 

dw  du 
{~dy  dz 

The  third  term  on  the  r.h.s  of  equation  (3.19)  is  known  as  solenoid  (or 
baroclinic)  term  of  the  circulation  equation  (3.13)  applied  to  an  infinitely  small 
horizontal  circuit. 

The  last  term  on  the  r.h.s  of  equation  3.19  is  the  frictional  term  ( or  also  beta 
term)  that  arises  due  to  turbulence  in  the  atmosphere. 

Calculation  from  synoptic  data  indicates  that  the  contribution  of  divergence 
term  [(f+  C )D]  is  more  than  the  other  terms  in  large  scale  weather  systems. 

So  if  little  divergence  exists  and  other  terms  are  negligible 

(3-20) 

Which  means  (/  + Q is  constant.  That  is  in  the  case  of  a ring  of  air  parcel 
revolving  round  the  earth,  the  absolute  vorticity  is  constant.  Which  means  as  ‘f 
increases,  the  relative  vorticity  (£)  should  decrease  and  vice  versa  in  order  that 
(/  + 4”)  to  be  constant.  So  if  air  moves  pole  ward,  ‘f  increases  and  so  Q decreases. 
Which  means  pole  ward  moving  air  must  lose  cyclonic  curvature  and  tend  toward 
anti-cyclonic  motion.  The  converse  is  true  for  the  air  moving  equator  ward.  Thus 
the  trajectory  of  an  air  parcel  in  the  meridional  direction  can  be  drawn  using  the 
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constant  absolute  vorticity  (CAV)  values  plotting  on  a map  as  shown  in  Fig.  11.  So 
for  weather  forecasting  CAV  trajectories  are  useful. 


Fig.ll.  CAV  trajectories  in  the  Northern  hemisphere 


Further  if  the  depth  of  the  column  of  the  air  parcel  ‘D’  is  considered,  the 
vorticity  theorem  will  be  as  given  in  equation  3.20. Where  ^ jis  called  the 
potential  vorticity  which  is  to  be  conserved  instead  of  simply  absolute  vorticity. 


Chapter  -4 

Quasi  Geostrophic  theory 

4.1.  POTENTIAL  VORTICITY  THEOREM: 


For  the  case  of  incompressible  baro tropic  fluid,  the  potential  vorticity 
equation  can  be  derived  using  the  continuity  equation  (2.3): 


& Sv  j . . . . ( du  dv\  dw 

1 1 = 0 which  means .’. 1 = . . 

dx  dy  dz ) ydx  dy ) dz 


..(3.23) 


The  vorticity  equation  (3.19)  can  be  written  with  the  help  of  equation  (3.23) 
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.(3.24) 


In  the  case  of  barotropic  condition,  the  thermal  wind  vanishes  and  so  the 
geostrophic  wind  is  independent  of  height.  Then  integrating  vertically  from  z=Di 
to  D2  equation  (3.24)  can  be  written  as 

+((+f)dz  = (f+()dw  on  integration  gives  rise  to 
dt 

=Cf 


r {C  + /) g\D2  - D, ] = {f  + C)g [WD2  - Wm ] (3.25) 


The  suffix  ‘g’  denotes  under  geostrophic  condition.  But  since  W = — we  can 

dt 


further  write 


Wd2  = and  Wdi  = ^1.  and  let  H =D2  -Di  is  the  depth 
dt  dt 


of  the  column  we  can  write  (3.25)  as 

d_ 
dt 


(c+f),[D2-DMf+c\[^-^]  - (/+<r)|(z>2-A) 


H.  (C  + ./  ) = l / +f ) d (H ) re  arranging  the  terms 
dt  dt 


( — - — )—  (C  + f)  = — — (h)  which  gives  rise  to 

C + f dtKh  ’ H dty  ’ 6 


d n 

— [b 
dt 


d 

or  — 
dt 


C + f 

H 

C + f 

H 


] = 0 


= 0 This  is  the  potential  vorticity  equation. 


If  the  flow  is  purely  horizontal  (w  = o),  equation  (3.24)  reduces  to 


- °- 

This  is  called  the  barotropic  vorticity  equation,  which  states  that  the  absolute 
vorticity  is  conserved  which  means  (C,  + /)  = constant. 
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4.2.  THE  BASIC  EQUATIONS  IN  DIFFERENT  VERTICAL  COORDINATES: 


Earlier  we  have  seen  the  equations  of  motion  in  Cartesian  and  spherical  polar 
coordinate  systems  taking  z as  vertical  coordinate.  But  modelers  find  it  convenient  to  use 
the  vertical  coordinate  as  pressure,  potential  temperature  (entropy)  or  the  vertical  stability 
parameter  (a)  and  derive  the  equations  of  momentum,  continuity,  vorticity  etc.  As  density 
does  not  appear  in  these  equations,  the  complicacy  of  its  variation  can  easily  be 
eliminated  by  expressing  these  equations  in  a modified  vertical  coordinate. 


4.3.  TRANSFORMATION  EQUATIONS:  ISOBARIC  COORDINATE  SYSTEM 

Since  pressure  and  height  are  directly  related  in  the  atmosphere,  it  is  common  to 
express  the  variables  in  terms  of  pressure  surfaces  rather  than  geometrical  heights  in 
meteorology.  So  it  is  essential  to  replace  the  vertical  coordinate  (z)  by  pressure  (p).  Once 
pressure  is  taken  as  vertical  coordinate  instead  of  z,  it  is  called  isobaric  coordinate 
system. 

When  ‘p’  is  taken  along  ‘z’  direction,  the  pressure  coordinate  is  not  normal  to  the 
isobaric  surfaces  but  is  exactly  vertical  upward  on  the  plumb  line. 

Consider  a cross  section  in  x-z  plane  and  let  some  arbitrary  property  ‘q’  is  constant 
along  the  small  section  of  the  surface  at  an  instant  of  time  as  shown  in  Fig.  12.  Let  the 
variable  is  £(x,z,t). 


Then  we  can  write 

(^')1_2=(^)1_3+(^)2_3 

for  small  increments  we  can  write 
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— — — = — — — + — — — — which  can  be  written  for  small  increments  as 
Ax  Ax  Az  Ax 


As  we  can  use  any  variable  in  place  of  C,  we  can  write  (1. 1)  as 


Suffix  ‘q’  implies  ‘holding  q as  constant’  in  the  operation. 

For  isobaric  coordinate  system  q = p(pressure).  A derivative  with  respect  to  z can  be 

written  as  JL  = ^ JL  ancj  if  we  make  the  hydrostatic  assumption  this 

dz  8z  dp 

becomes  — = -pg  — (4.3) 

dz  dp 


So  using  (4.3),  equation  (4.2)  may  be  written  as 


Similarly  we  can 

write  for  ‘y’ 

and  ‘t’  as 

(d\  _ 

{dyjz 

r a 
l^yJ; 

.*“(1), 

5 

dP 

(4.5) 

f5l  -1 
[dt  J,  1 

-1 

,dt)p 

*411 

d 

dp 

(4.6a) 

The  last  term  of  equation  (4.6a)  denotes  the  rate  at  which  the  height  of  an  isobaric 
surface  changes  with  time.  Equations  (4.4)  to  (4.6a)  denote  the  inputs  to  transform  the 
(x,y,z,t)  Cartesian  coordinate  system  to  (x,y,p,t)  isobaric  coordinate  system. 


4.3.1.  ADVANTAGES  AND  DISADVANTAGES  OF  ISOBARIC  COORDINATE 
SYSTEM: 


The  isobaric  coordinate  system  has  several  advantages  for  modeling  work.  Some  of 

them  are: 

a)  meteorological  data  is  normally  referred  to  isobaric  surfaces. 

b)  the  continuity  equation  in  this  system  is  reduced  and  in  simple  form 

c)  density  does  not  explicitly  appear  so  no  problem  of  integration  or  differentiation 

d)  sound  waves  are  completely  filtered. 

The  serious  disadvantage  of  this  system  is  found  in  modeling.  Modellers  use 
boundary  conditions.  Here  the  lower  (z0)  boundary  condition  is  assumed  coincident  with 
the  surface  pressure  (p0).  we  know  surface  pressure  is  not  constant  on  the  level  or 
geopotential  surface.  Even  when  the  ground  is  level  also  surface  pressure  changes  in  the 
horizontal. 
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4.3.2.SIGMA  COORDINATE  SYSTEM: 


To  avoid  the  serious  disadvantage  of  the  isobaric  coordinate  system,  cr 
coordinate  system  has  been  introduced  in  modeling  work.  The  disadvantage  of  isobaric  - 
coordinate  system  is  the  coincidence  of  surface  pressure  at  the  lower  boundary  with  the 
ground  (z0)  which  is  not  true. 

In  the  sigma  ( cr ) system,  the  vertical  coordinate  is  the  pressure  normalized  with  the 
surface  pressure  (ps). 

:.a  = ^~  (4.6b) 

Ps 

This  is  called  the  static  stability  parameter.  Thus  cr  here  is  a non  dimensional 
independent  vertical  coordinate  which  decreases  upward  from  a value  cr  = 1 at  the 
ground  to  cr  = 0 at  the  top  of  the  atmosphere.  Here  in  cr  coordinates  the  lower  boundary 
condition  will  always  apply  exactly  at  cr  =1.  Further  the  vertical  sigma  velocity  is  given 
as 

da 

a = — 
dt 

This  will  always  be  zero  at  the  ground  even  in  the  presence  of  sloping  terrain. 

So  the  lower  boundary  condition  in  the  cr  system  is  cr  = 0 at  cr  = 1.  Thus  this  cr 
coordinate  system  is  particularly  useful  in  the  regions  of  strong  topographic  variations. 


4.4.  OMEGA  EQUATION: 


We  know  the  differential  variable  in  x,y,z,t  coordinate  system  as 

d d 8 8 8 

— = — + u — + v — + w — (4.7) 

dt  dt  8x  dy  dz 

For  converting  this  to  isobaric  coordinate  system,  apply  equation  4.6  to  4.7  then 


(d 


t:  = hr  +«hr 


a 


v'-l  +P8\jt)p-dp’ 


..(4.8) 


[v  comparing  equations  4.7  and  4.6,  w—  = pg\—  j —and  we  know  — = w\ 

dz  \8t  J dp  dt 


Further  equation  (4.8)  can  be  written  as 
d_ 
dt  " 


1 f 8 1 

dpi 

' d_A 

Uu 

J,  “U, 

l+vUJ, 

dz\ 

,dtjl 

dp 


_dp 


We  can  write  equation  (4.9)  taking  <y(omega)=  (4.9) 

dt 


as 

d f 0 ^ f 0 ^ dp  0 

— = — +u\ — + v — — 

dt  \dtjp  \dxjp  \8y J p dt  dp 


.(4.10) 
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Here  the  variable  &>(omega)=  - —is  equivalent  to  vertical  motion  which  is  like 
dt 

— = win  x,y,z,t  coordinate  system.  The  negative  sign  implies  that  the  rate  of  pressure 
dt 

decreases  with  height.  So  (4.10)  becomes 


d 

). — . 
dp 


..(4.11) 


This  is  called  the  omega  equation.  Thus  the  individual  derivative  in  isobaric 
(x,y,p,t)  coordinate  system  takes  the  same  form  as  in  Cartesian(x,y,z,t)  system  except  that 

co.— replaces  w — . So  the  variable  co  is  equivalent  to  vertical  velocity  in  isobaric 

dp  dz 

coordinate  system.  When  measured  in  microbars  per  second,  co  is  numerically  equal  to  w 
in  cm  s'1  in  the  lower  atmosphere  with  opposite  sign. 


4.5.  MOMENTUM  EQUATION  IN  ISOBARIC  COORDINATE  SYSTEM: 


We  know  in  Cartesian  (x,y,z,t)  coordinate  system  momentum  equations  as 

(4.12) 


du  du  du  du  du  dp  , 

— = h u hv h w — = -a 1-  jv 

dt  dt  dx  dy  dz  dx 


dv  dv  dv  dv  dv  dp  r 

— = — + u — + v — + w — = — a — ~ju 
dt  dt  dx  dy  dz  dy 

In  a similar  way  the  momentum  equations  in  isobaric  (x,y,p,t)  coordinate  system 
can  be  written  from  equation  (4.10)  by  incorporating  u and  v as 


(du 
u\  — 

lac 


(du\ 
f v — 


f _( ( dv ) [ dv  | 

\dtjp  \dt)p  \dx)  p 


du 

dz 

,+(°'dp ~ 

8 dx 

dv 

dz 

+ 0}'dp ~ 

8 dy 

+ fv  (4.12a) 


..(4.12b) 


Please  note  adp  = gdz  is  used  here  to  the  first  term  of  the  right  hand  side.  This  kind 
of  momentum  equations  are  used  by  modelers  because  the  density  term  does  not  appear 
explicitly. 


4.7.  THE  VORTICITY  EQUATION  IN  ISOBARIC  COORDINATE  SYSTEM: 


This  form  of  vorticity  equation  can  be  obtained  by  differentiating  momentum 
equation  in  isobaric  coordinates  (4.12a)  with  respect  to(  w.r.t)  ‘y’  and  the  equation 
(4.12b)  with  respect  to  ‘x’  and  subtract  the  former  from  the  latter  holding  ‘p’  as  constant. 
Then  we  get 


.(4.12a) 
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.(4.12b) 


Y)  +/Y|  +/* 

dt)p  ydx)p 

Differentiating  (4.12a)  w.r.t  ‘y’  we  get 

& 

dydt 


dv  f dz)  r 

(o.-  = -g\  — -fu 
dp  ^qy  J 


8u  du  ( d2u  | du  dv  d2u  du  dco  _ d2z 
dx  dy  dy  dy  dydp  dp  dy  ^ dxdy 

(4.19) 


Differentiating  (4.12b)  w.r.t  ‘x’  we  get 


J [axdy) 

entiating  (4. 1 

d2 3v  ^ ( a2v^l  dv  du  ( d2v  ^ dv  dv  d2v  dv  dco  d2z  ..  du 

dxdt)  [ dx  J dx  dx  \ dxdy  j dy  dx  dxdp  dp  dx  dxdy  dx 

Subtracting  (4. 19)  from  (4.20)  and  taking  ( — - — I = ^ = k .( V .V)  (we  get 

dy  1 

J V^J  ^ayj  dp  dy) 


(4.20) 


^ ^ , dC  J du  dv\  dv  dco  du  dco  J du  dv 

l)t}  + U{lk}  + V{~ty  ] + (0'~dp+\~d^  + dyj  + ty'lh-ty'lty-~f{lk  + dy 


(d£  ) 

{ dt  ) 


f 


4 


-v(3- 


dad | Sv  ( dco  j du 
,dx)pdp'  { dy)pdp ' 
(4.21) 


This  is  the  vorticity  equation  in  isobaric  coordinate  system.  The  first  term  on  the 
r.h.s  is  divergence  term,  the  second  term  is  /?  plane  term  and  the  last  two  terms 
considering  as  the  third  term  is  tripping  term.  Note  that  the  solenoidal  term  is  missing 
here.  This  is  because  the  isobaric  system  doesn’t  depend  on  the  density  varaiation. 

This  can  be  further  written  taking  — = 0 as  the  coriolis  parameter  (f)  at  any  latitude  is 

dx 

constant  and  — = J3 , variation  of  coriolis  parameter  with  latitude. 

dy 

5£  = -V.V(f  + /)-*.«£ _(f  + f)vy  - Where  V = iu  +jv (4.22) 

dt  dp  dp  dx  dp  dy 


4.8.  QUASI-GEOSTROPHIC  VORTICITY  EQUATION: 

Referring  to  the  vorticity  equation  in  isobaric  coordinate  system  in  equation  (4.22),  the 
terms  reading  from  left  to  right  are  as  follows: 

1 . The  local  rate  of  change  of  relative  vorticity 

2.  the  horizontal  advection  of  absolute  vorticity 

3.  the  vertical  advection  of  relative  vorticity 


dx 


df  df 

u--v— 
dx  dy 
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4.  the  divergence  term 

5.  the  twisting  or  tilting  or  tripping  term 

As  shown  by  the  scale  analysis  earlier  we  simplify  the  vorticity  equation  for  midlatitude 
synoptic  scale  motions  by 

1 . neglecting  the  vertical  advection  and  twisting  terms 

2.  neglecting  C,  compared  to  f in  the  divergence  terms 

3.  approximating  the  horizontal  velocity  by  the  geostrophic  wind  in  the  advecting 
term  and 

4.  replacing  the  relative  vorticity  by  its  geostrophic  value  and 

5.  expressing  coriolis  parameter  in  Taylor’s  series  with  initial  value  f0  at  latitude  tp0 
such  that 


f = f0  + Py  + higher  order  terms  (4.23) 

where  P = — and  y = 0 at  cp0-  Neglect  the  higher  order  terms. 


P plane  approximation: 

Let  L designate  the  latitudinal  scale  of  motion,  then  the  ratio  of  the  first  two  terms  of  ‘f 
has  the  order  of  magnitude  as 

L 

J3L  _ dy  _ d(2Qsm<f>0)L  _ (2Qcos^0)L 
fo  fo  2Qsin  fa.a  2Qsin  </>0.a 
Thus  when  the  latitudinal  scale  (L)  of  motions  is  small  compared  to  the  radius  of  the 
earth  (a),  (then  L<a)  we  can  let  the  Coriolis  parameter  f0  has  a constant  value  except 

where  it  appears  differentiated  in  the  advection  term.  In  such  case  — = P is  assumed  to 

dy 

be  constant.  This  approximation  is  usually  referred  to  as  the  beta-palne  approximation. 
Applying  the  above  assumptions  and  approximations,  equation  (4.22)  becomes: 

^ = -Vg  .V(^  + /)  - (/0  )V.V  (4.24) 


Where  £ = V 9 and  V =kx^ 
g fo  S fo 

This  is  called  quasi-geostrophic  vorticity  equation  . 

It  is  important  to  note  from  equation  (4.24),  in  the  divergence  term  horizontal  wind  is  not 
replaced  by  its  geostrophic  value.  This  is  because  the  small  departures  of  the  horizontal 
wind  from  geostrophy  accounts  for  the  divergence. 

We  know  the  continuity  equation  from  equation  (4.17) 
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Substituting  this  equation  (4.25)  in  (4.24)  we  get 

^ = -V„.V(fI+/)  + (/0)^  (4.26) 

dt  dp 

This  is  the  alternative  form  of  quasi  geostrophic  vorticity  equation.  The  advantage  of  this 
equation  is  a>  estimates  are  better  than  the  estimates  based  on  continuity  equation  (4.17) 
because  both  the  local  change  of  geostrophic  vorticity  and  the  advection  of  vorticity  can 
be  estimated  accurately  as  they  are  related  to  geopotential.  This  equation  (4.26)  forms  as 
a closed  prediction  equation  in  co  . 


28 


